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Abstract 

The double J /if) production in pp collision is revisited. It is found that the J/if> 
pair production rate at leading order with px > 4 GeV in conventional scheme 
is larger than the rate derived from color-octet mechanism by a factor of five. 
As the double J /if) production is shown to be attainable with data yet collected 
at the Tevatron detectors, experimental study on it would be helpful to clarify 
the validity of the conventional color-singlet description for heavy quarkonium 
production. 
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Quarkonium production and decays have long been taken as an ideal means to investi- 
gate the nature of QCD and other new phenomena. Hence, to establish a proper theory 
which can precisely describe heavy quarkonium production and decays is very necessary. 
A novel effective theory, non-relativistic QCD(NRQCD) is possibly one theory to this 
aim which is formulated from the first principles. However, to make a precise prediction for 
quarkonium production, with only the NRQCD is not enough at least for now, since the 
magnitude of the non-perturbative parameters in the theory are still unknown. Therefore, 
generally speaking theoretical calculation results may depend on the values of the input 
non-perturbative parameters, although NRQCD can give out the relative weights of these 
parameters in order of v 2 from its "velocity scaling rules". Up to now, on one hand the 
color-octet mechanism still stands as the most plausible proposal in explaining the large 
transverse momentum production "anomaly" ||; on the other hand, it encounters 

some difficulties in confronting other phenomena ||. Especially, a recent research || shows 
that previous studies of the color-octet contributions to high-p T J/ipijp') production were 
overestimated, although not much, which may somehow bring new light to the understand- 
ing of heavy quarkonium, especially for charmonium, production mechanism. All in all, to 
what degree the color-octet mechanism plays the role in quarkonium production is still not 
clear and an interesting question. 

During past decade, lots of efforts have been made in proposing approaches to find definite 
color-octet signals. Unfortunately, the nature appears to be more elusive than expected and 
experiments still give no conclusive answer to this demand yet. One interesting point is 
that once there is a proposal highlighting a octet process, not very long people find another 
competing process in singlet, which brings difficulties to experimenters in distinguishing 
between color-octet and singlet. For instance, in electron-position scattering Ref. to Refs. 
0, and Ref. || to Ref. || in direct photon-photon collision. 

In explaining the high-pr ^ surplus production discovered by CDF group [§] at the 
Fermilab Tevatron, the color-octet scheme tells us that the dominant source of charmonium 
production at large transverse momentum comes from the produced hard gluon followed 
by its fragmentation into an intermediate color-octet cc state, which eventually evolves into 
quarkonium non-perturbatively. Applying this idea to the double-quarkonium production, 
Barger et al. |IJ studied the quarkonium pair production via the double-gluon fragmentation 
and found that this scheme for double J/ip production gives a result which is measurable 
with the accumulated data at the Fermilab Tevatron detectors. And hence, claimed the 
measurement of this process in experiment would give a test on the color-octet mechanism. 
At sufficient high transverse momentum this treatment might be true, however, at moderate 
and relatively low px one should take great care. In the present work we show that it is 
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really the case for double- J/ip production at even not very low transverse momentum. 

In color-singlet model, the partonic subprocess starts at order a 4 . The lowest order 
processes include g + g —> J/tp + J/tp and q + q — > J/tp + J/tp. Since the latter, the quark 
annihilation process, reasonably gives less contribution at Tevatron energy, in this brief 
report, we restrict our calculation only to the gluon-gluon fusion one, as shown in figure 
1. This parton process is similar to the case in photon-photon scattering discussed in Ref. 
||, however, the extension to the present discussion is not trivial. The QCD non-Abelian 
nature makes more possible channels involve in, i.e., the third topological group shown in 
the Figure. 

The differential cross section for quarkonium pair hadroproducion is given by 
^-{pp -> U/tp + X) = J2j dx a d yi f a/p (x a )f b/p (x b ) 2 ^ T _ X f T b eyi ^(a + b - 2J/^ + X) , (1) 



a.b 



where f a / p and fyp denote the parton densities; t, as well as s and u appearing later, is 
Mandelstam variable of the parton level; 2/1(1/2) is of the rapidity of produced J/tp; x T = ^jX- 
with rriT = \fm? + p\. Here, the capital vS denotes the incident beams total energy, and 
taken to be 1.8 TeV of the Tevatron in our subsequent numerical calculations. 

The parton scattering differential cross section of gluon-gluon to J / tp pair can be calcu- 
lated by the standard way straightforwardly. It is 

da _ 16 a* n \R(0)\ 4 

dt 81m 2 s 8 (m 2 — t) 4 (m 2 — u) 



[2680m 24 - 14984m 22 t + 31406m 20 t 2 - 31824m 18 t 3 



+ 17668m 16 t 4 - 7172m 14 t 5 + 2956m 12 t 6 - 794m 10 t 7 + 47m 8 t 8 + 20m 6 t 9 + m 4 t 10 

- 14984m 22 w + 89948m 20 tw - 202576m 18 t 2 w + 228560m 16 t 3 w - 153360m 14 t 4 w 

+ 76406m 12 t 5 w - 30782m 10 t 6 w + 7642m - 822m 6 t 8 w - 66m 4 t\ + 31406m 2 V 

- 202576m 18 tu 2 + 470856m 16 tV - 536476m 14 tV + 361624m 12 tV - 182454m 10 tV 
+ 73146m 8 tV - 17902m 6 tV + 2469m 4 t 8 w 2 + 36m 2 tV - 31824m 18 w 3 + 228560m 16 tu 3 

- 536476m 14 tV + 571900m 12 t 3 w 3 - 335186m 10 tV + 150334m 8 tV - 58126m 6 t 6 u 3 
+ 12874m 4 tV - 2344m 2 tV + 17668m 16 w 4 - 153360m 14 tw 4 + 361624m 12 tV 

- 335186m 10 tV + 132502m 8 tV - 35306m 6 tV + 11928m W - 148m 2 tV + 698*V 

- 7172m 14 w 5 + 76406m 12 tw 5 - 182454m 10 tV + 150334m 8 tV - 35306m 6 tV + 1164m 4 tV 

- 1576m 2 tV - 1816tV + 2956m 12 u 6 - 30782m 10 tu 6 + 73146m 8 tV - 58126m 6 tV 
+ 11928m 4 t 4 u 6 - 1576m 2 t 5 w 6 + 2748t 6 u 6 - 794m 1 V + 7642m 8 tu 7 - 17902m 6 tV 

+ 12874m 4 t 3 w 7 - 148m 2 tV - 1816*V + 47mV - 822m 6 tw 8 + 2469m 4 tV 

- 2344m 2 tV + 698tV + 20m V - 66m 4 tw 9 + 36m 2 t 2 u 9 + m 4 u w } , (2) 
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where m is the mass of charmonium, the J/ip; \R(0)\ is the magnitude of its radial wavefunc- 
tion at origin. In obtaining the above analytical expression, we start from general Feynman 
rules and project the Charm-anti-Charm pair into the S-wave vector charmonium state in 
color-singlet. To manipulate the trace and matrix-element-square of those tens of diagrams, 
the computer algebra system MATHEMATICA is employed with the help of the package 
FEYNCALC (TJ. 



The values of input parameters used in our numerical calculations are: 

m c = 1.5 GeV, |R(0)| 2 = 0.8 GeV 3 . (3) 

Here, the non-relativistic relation m = 2m c is adopted. The typical scale is set to be at mr, 
and hence the strong coupling is running with transverse momentum. 

With the above preparation, direct J/0 pair production cross section can be imme- 
diately obtained. Applying the pseudorapidity cuts on both produced charmonia, i.e. 
|?7(^i)|, 1 77 ('02) I < 0-6, the integrated cross section a{pp — > Vy^-Vy^-) for pr(0) > 4 GeV is 
0.70 pb. Here, the notation "0^+^- means that the branching ratio of B(ip — > = 0.0588 

of practical measuring mode to reconstruct the charmonium state is included. For integrated 
luminosity of 100 pb^ 1 achieved in the past run of the Tevatron, there should be about sev- 
enty J I ip pair events coming from the conventional scheme, which is about five times larger 
than what obtained in Ref. | ID| . In our numerical calculation, the parton distribution of 



CTEQ5L |T!| is used, and both renormalization scale and factorization scale are evolved to 
the same point my. In the numerical calculation, the integration limits of x a in Eq. ([[]) are 
truncated in accordance with the pseudorapidity cuts in both produced charmonia. 

The double- J/tp production cross section as function of transverse momentum p? is shown 
in Figure 2 in solid line. For comparison, the result from color-octet process discussed in Ref. 
|j~0|| is presented in the same figure in dashed line. From the figure we see the conventional 
production scheme dominates over the octet one in relatively low-p^ region, p? < 7 GeV, 
where most of the charmonium pairs are produced. In principle, the color-octet mechanism, 
at least the importance of it, can be undisputedly checked in large transverse momentum 
area, however, unfortunately there is no enough available data for such purpose at the 
moment. 

In above analysis, we take charmonium state, the J/ip, as an object. Nevertheless, the 
results can be readily applied to some other quarkonium states, both charmonium and bot- 
tomonium. For example the ■0' pair production rates in color-singlet model can be obtained 
by multiplying the constant 



li^O)! 4 B 2 {J/i) -> fi+fi-) B\Jji> -> fi+fi-) 



(4) 
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and the J/ip results. 

In conclusion, we have calculated the J / ip pair production rate at the Fermilab Tevatron 
with the conventional heavy quarkonium production treatment, the color-singlet model. The 
calculation is carried out at the leading order in strong coupling constant. Whereas, the 
higher order corrections can be taken to be small due to relatively high interacting scale we 
have considered; the relativistic correction, which is not always too small to be negligible in 
dealing with charmonium system, is also double suppressed since we are considering the pair 
production; and for the same reason, the higher excited states feeddown is well suppressed as 
well. From our analysis, the accumulated data at the Fermilab Tevatron detectors can give 
about one hundred of J/tp pair events with transverse momentum p T > 4 GeV, among them 
about eighty percent is produced via the conventional quarkonium production mechanism. 
This finding tells us that analyzing the data for this process is not a suitable game in 
chasing the color-octet signatures, but rather provides an ideal means to precisely check the 
conventional description for heavy quarkonium production. 
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FIGURES 




FIG. 2. The differential cross-section of J ftp pair production versus pt at the Tevatron. Solid 
line comes from the Color-Singlet calculation of this paper; the dashed line from the Color-Octet 



calculation read from Ref. [10|. 
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